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Abstract 

We consider the backreaction problem of a quantized minimally coupled massless scalar field in 
cosmology. The adiabatically regularized stress-energy tensor in a general Friedmann-Robertson- 
Walker background is approximately evaluated by using the fact that subhorizon modes evolve 
adiabatically and superhorizon modes are frozen. The vacuum energy density is verified to obey 
a new first order differential equation depending on a dimensionless parameter of order unity, 
which calibrates subhorizon/superhorizon division. We check the validity of the approximation by 
calculating the corresponding vacuum energy densities in fixed backgrounds, which are shown to 
agree with the known results in de Sitter space and space-times undergoing power law expansions. 
We then apply our findings to slow-roll inflationary models. Although backreaction effects are found 
to be negligible during the near exponential expansion, the vacuum energy density generated during 
this period might be important at later stages since it decreases slower than radiation or dust. 



* ali.kaya@boun.edu.tr 

^ merve.tarman@boun.edu.tr 



1 



I. INTRODUCTION 



Quantum field theory in curved space is a well established subject of study (see e.g. [1- 
3]). For free fields in a general background, difficulties related to the absence of Poincare 
symmetry, non-uniqueness of vacuum and lack of particle states were all solved and a physi- 
cally viable picture has been obtained. Suitable regularization techniques, like adiabatic [5] 
(see also [6]) and point-splitting [4], were developed to obtain finite and physically relevant 
expressions for the stress-energy-momentum tensor vacuum expectation values. Adiabatic 
and point-splitting regularizations were shown to be equivalent in cosmological spacetimes 
[7, 8]. 

Although quantization in a fixed background is well understood, the backrcaction of a 
quantum field on the geometry is not fully determined. For example, the cosmic evolution 
of the scale factor of a universe dominated by the stress-energy tensor of a quantum field is 
not known. The backreaction problem is difficult to tackle (see e.g. [9-14] for some attempts 
relevant for cosmology) and it is expected to shed further light on semiclassical quantum 
gravity. 

In this paper, we try to determine the cosmological backreaction of a minimally coupled 
quantized real massless scalar field. As we will see, after adiabatic regularization the problem 
reduces to an integro-differential equation system involving the scale factor of the universe 
a{t) and the mode function of the field fJ,k{t)- In an earlier work [15], we simplified this 
equation system by using adiabatic subtraction terms up to second order and approximately 
evaluating the momentum integrals from the behavior of the subhorizon and the superhorizon 
mode functions. In that way we managed to evaluate the mode function momentum integrals 
in terms of the vacuum energy density and obtained analytical solutions. 

We would like to improve the method employed in [15] in various ways to have a better 
approximation. First, we utilize the full adiabatic regularization, i.e. we include adiabatic 
order four subtraction terms to remove all the divergences in the quantum stress-energy 
tensor. Second, we introduce a dimensionless parameter a of order unity characterizing 
subhrorizon/super horizon division of modes. Finally, we also keep finite adiabatic subtrac- 
tion terms, which vanish when one naively takes the massless limit before evaluating the 
momentum integrals (for example these terms must be kept to obtain the trace anomaly 
using adiabatic regularization [16]). As a result, the final equation system involving the 
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scale factor and the vacuum energy density turns out to be more complicated than the one 
found in [15] and it is difficult to obtain analytical solutions. However, one can still obtain 
quantitative information on how vacuum energy density affects the cosmic evolution. 

To check the validity of our approximation we determine the vacuum energy density in 
a fixed background to compare our findings with the known results in de Sitter space and 
space-times undergoing power law expansions. We find complete agreement in both cases 
by adjusting the above mentioned dimensionless parameter a. We also apply our formulas 
to slow-roll inflationary models and try to pin down possible backreaction effects. We show 
that the vacuum energy density becomes negligible during the near exponential expansion, 
but it may have a cosmological imprint at later stages since it decreases more slower than 
radiation or dust. 



II. QUANTIZATION AND REGULARIZATION 

We consider a free massless real scalar field obeying 

W = 0. (1) 

The scalar is propagating in a Priedmann-Robertson- Walker (FRW) space-time which has 
the metric 

ds^ = -dt^ + a{tf{dx^ + dy'^ + dz^) 

^a{rif{-drj^ + dx^ + dy'^ + dz^), (2) 

where t and r] denote the proper and conformal cosmic times, respectively. We denote the 
corresponding Hubble parameters as 

h=-, H=-, 3 
a a 

where the prime and the dot represent derivatives with respect to rj and t, respectively. 
Defining a new field by 

II = acf) (4) 

and using the standard canonical quantization procedure with the action S — J \f—g (V0)^, 
one can decompose the field operator as 



-* J. 

where k is the comoving momentum variable, ar and al denote the time-independent lad- 



der operators obeying [a^, at ] = 5{k — k'). Furthermore, the mode functions satisfy the 
Wronskian condition 

i^kiJk - l^ll^'k = h (6) 

together with 



k' 



a 

a 



Hk = 0. 



The ground state |0 > at time 770 can be defined by imposing 



(7) 



(8) 



The quantum system is uniquely specified if one additionally gives the initial values iJ>k{vo) 
and fi'kivo) obeying the Wronskian condition (6). 

Using the definition of the stress-energy-momcntum tensor Tf^^ = V ^(pVpCt) — ^gij.i^(V4>y, 
one can calculate the following vacuum-expectation values 
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These expressions contain quartic, quadratic and logarithmic divergences. Adiabatic reg- 
ularization [5] offers a physically viable way of removing these divergences using suitably 
defined subtraction terms. To determine the adiabatic subtraction terms, one defines the 
adiabatic mode function n'jf in terms of a new variable Qk in the following WKB form 



Using (7), Qfe can be seen to obey 



a ^4Q| 



2VLk 



(10) 



(11) 



Eq. (11) can be solved itcrativcly, where the number of time derivatives acting on is used 
as an expansion parameter. The zcroth order adiabatic solution is fi^ = k and using it in the 
right hand side of (11) the second order adiabatic solution can be found as 0^ = /c 1 — . 
To regularize (9), i.e. to remove quartic, quadratic and logarithmic divergences, one should 
determine Jljt up to adiabatic order four terms and subtract the corresponding stress-energy 
tensor expressions obtained from the adiabatic mode function (10). After a straightforward 



calculation, the regularized finite expressions for the vacuum expectation values can be found 
as 



Pv 
Pv 



47r2a4 Jo 



1 



477^0^ Jo 



















a' 









24A;3 



k _h'^ h' 



{3h^ + h'^ - 2hh" - Uh^h' + 2h"') 



dk, (12) 



dk. (13) 



One may check that pv and Py obey the conservation equation py + 3h{pv + Pv) — 
provided (7) holds. 

There is an important subtlety in using adiabatic regularization for massless fields. If 
one naively performs adiabatic regularization for a conformally coupled massless scalar field, 
the regularized stress-energy tensor vanishes identically, i.e. it is not possible to obtain the 
conformal anomaly. The correct procedure is to start with a massive field and then send the 
mass parameter to zero by carefully keeping the terms that survives the limit. Only in this 
way one can get the conformal anomaly using adiabatic regularization [16] . In our case, the 
same procedure should be used, this time for a minimally coupled field. These terms can 
be determined by using the formulas given in [17]. Similarly, in [18] we also find out the 
adiabatic subtraction terms for a massive field up to sixth order. From our calculation in 
[17] we note the fourth order adiabatic subtraction terms for energy density as 

roo rlh r 

/ 56m^a;'a^a''a" - 105m«a^a'^ + UAm^Ma'''a" + 16m'a;«a=^(a"' - 2a'a"') 

Jo 5127r^a'co'^^ l 

+6Au^a'^a" - 22Am'^uj^a'a'^ + AmWa^'ia"^ - 2a' a'") + Uuj^'^aa'^ + IQm^uj^aa'* + U^lo^^o" 
+16m^a;^a^a'2 + lQuj^a{a"^ - 2a' a'") + 64m'a;«a=^a'' + l2Am''Ma'^ + QGrn^c^^a^a'^a" 



where oj = \/k'^ + m?a^. By scaling k — )> mk the integrals of the terms in the square brackets 
with m^, m^cu^, rn^u'^ and vn^uf' can be seen to be m-indcpcndent and finite (note that the 
integral of lo^ term is also mass independent but it diverges). Performing the integrals one 
can find the finite adiabatic subtraction terms for energy density 

PF = ^ \92h^ + GOh^h' + llh'^ - 22hh"\ . (14) 

The corresponding expression for pressure can be determined from the stress-energy conser- 



vation p'p + 3h{pF + Pp) = 0, which gives 

1 



28 807r2a4 
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(15) 



These are the extra terms which must be subtracted from (12) and (13) 



Pv^Pv- Pf- 



Pv ^ Pv - Pf 



(16) 
(17) 



One can use pv and Py as sources in the Einstein equations: 



(18) 




(19) 



The evolution is fully specified once the initial conditions for the mode function /Xfe and the 
scale factor a{t) are given. This is an integro-differential equation system: the integrals of 
the mode functions determine the stress-energy tensor as in (12) and (13), and these are 
used in Einstein equations to determine the dynamics of the scale factor. We have also (7), 
which fixes the evolution of the mode functions. 

Although (12) and (13) are guaranteed to be ultraviolet finite for suitably chosen mode 
function initial conditions, there are potentially dangerous infrared divergences for massless 
fields [19]. For example, when a = {tiq/viY the mode function for the Bunch-Davies vacuum 



is given by Hk = t] ^ k / 2 hn{kri) , where hn is the spherical Hankel function of first kind. 
From the behavior of hn{kri) as A; ^ 0, one can see that the integrals of the mode functions 
in (12) and (13) diverge with a power near k — Similarly, the integrals of the fourth order 
adiabatic subtraction terms in (12) and (13) (these are the terms with in the square 
brakets) diverge logarithmically near k — Therefore, one should introduce an infrared 
cutoff to make sense of pv and Py- This is similar to the point splitting regularization, 
which also requires an infrared cutoff for massless fields. 

III. AN APPROXIMATION 

In this section we would like to utilize an approximation to simplify the integrals in (12) 
and (13). First, we assume that the field is placed in Bunch-Davies vacuum or Tth order 
adiabatic vacuum with I > 4 such that as A; — > oo one discovers usual mode functions of the 
Minkowski space. ^ The differential equation determining the evolution of the mode functions 

^ One should note that in reahstic situations there is a certain uncertainty in specifying the vacuum [20]. 
See also [21] for possible cosmological implications of different initial states in inflation. 
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(7) can be solved in two different limits: Hk — e^'^'^/V2k ii k ^ h and (fjik/o)' ^ ii k <^ h, 
respectively, corresponding to the subhorizon and the superhorizon regimes. 

To simplify (12) and (13) we use the fact that subhorizon modes evolve adiabatically and 
superhorizon modes freeze-out. Therefore, in the chosen vacuum the difference between the 
actual mode function /x/,. and the adiabatic mode function /x^'^ in (10) should be negligible 
for subhorizon modes. Aditionally, we take superhorizon modes to obey exactly {^k/a)' — 
0. To have a better approximation, we introduce a real parameter a, which quantifies 
subhorizon/super horizon border. Namely for k > ah we ignore the difference between /ik 
and iJ.l''' in the integrals in (12) and (13), and for k < ah we take (/ifc/a)' = 0. The 
need for such a parameter can be justified as follows. From the differential equation (7), 



sub horizon/super horizon regimes are determined by ya"/a, which may be different than h. 
Moreover, there appears errors in the integrals for using //^"^ instead of fik for subhorizon 
modes and imposing (iik/a)' — exactly for superhorizon modes. The a-parameter is 
introduced to compensate these errors and it is expected to be of order unity.^ 

One can evaluate (12) and (13) approximately as follows. First, the integrals are negligible 
in the [ah, oo] range by the above argument. To calculate the integrals in [0, ah] interval we 
define the following variables 
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Since we approximate superhorizon modes to obey^ il^k/o)' = 0, these variables satisfy 

Ps = -\ps. 

Performing the integrals (12) and (13) one finds 



Pv ^ Ps- 
Pv=Ps- 



A^:'^a'^ 
1 

4^ 



{a^ + a^) 
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{3h^ + h'^ - 2hh" 



- PF: 



(20) 
(21) 

(22) 

(23) 



'-h' - —h'h' 
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1 ^ / ah^ 

+24 ""l:^^ 



{3h^ + h'^ - 2hh" - 12h^h' + 2h"') 



-Pf 



(24) 



^ Our a parameter is similar to the e parameter of Starobinsky introduced in [22]. 

^ In the superhorizon regime, there are two solutions to (7): /J.kiv) — ci{k)a{r}) + C2(fc)a(r/) dr]' /a?{r]'), 
where (6) implies C1C2 — c*C2 = i. The second solution is the decaying mode and it becomes negligible in 
time. Indeed, {nk/a)' — 02/0^ and thus first terms in (20) and (21) decrease more, which justifies 
the equation of state (22). _ 



where a comoving infrared cutoff ho is introduced to make sense of the logarithmic integrals 
(the cutoff must be comoving in order to preserve stress-energy conservation). Infrared 
divergences, which might be related to the behavior of the mode functions, are absorbed in 

the definitions of ps and P5. Thus, the vacuum expectation value of the stress-energy tensor 
is determined up to an unknown function ps, which is actually fixed by the initially chosen 
vacuum state. 

The combination pv + SPy, which appears in the right hand side of the acceleration 
equation (19), becomes 

1 



Pv + 3Py = 



2 2 

-{Pf + 3Pf). (25) 



+ 1 In ( ^ ) {Qh' + 2/i'2 - 4hh" - 12h^h' + 2h"') 
8 \ho J 



Thus, Ps dependence disappears and (19) can be studied to determine the evolution of the 
scale factor 0(77). 

Note that pv and Py must obey stress-energy conservation. This is guaranteed for the 
exact expressions (12) and (13), and we can now impose it to determine the evolution of the 
vacuum energy density. Plugging (23) to the right hand side of p'y — —3h{pv + Py), using 
Ps — —ps/^ and further solving for ps in terms of pv from (23) one finds^ 



p'y + 2hp 



1 



V 



47r2 



(q;4 + a^) 



h' - ^h'h' 



2 



2 2 

+ h{pF + 3PF). (26) 



+ i In ( ^ ) {Qh' + 2hh'^ - Ah^h" + 2hh"' - 12h^h') 
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This is our final expression for the vacuum energy density. The backreaction problem is 
expressed in terms of two variables, namely the scale factor of the universe 0(77) and the 
vacuum energy density pv, which obey the Priedmann equation (18) and the conservation 
equation (26). As usual the acceleration equation (19) is satisfied identically once (18) and 
(26) are hold. 

Since pv obeys a first order differential equation, its initial value must be specified, which 
is fixed by the initial vacuum chosen. Introducing the homogeneous and particular solutions 
Pv — Ph + Pp, (26) implies that the homogeneous piece evolves as pn — Co/a?, where the 
constant Co can be related to the initial value of pv- The decrease of the homogeneous 



* Recall that pp and Pp are given in (14) and (15). These also obey stress-energy conservation, which we 
use in (26) to substitute for p'p. 
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piece is equivalent to a perfect fluid with an equation of state parameter ui — —1/3, which 
is slower than radiation and dust. 

In general, a parameter is expected to depend on time but its time dependence must 
be weak since a is of order unity. In any case, it is not possible to integrate (26) to fix 
pv exactly. However, one can estimate the magnitude of the vacuum energy density and 
determine how it evolves. 

The comoving infrared cutoff ho must also be chosen appropriately. One way of dealing 
with the infrared problem is to assume that the vacuum has no superhorizon modes initially 
(see e.g. [23]). Thus can be identified as the initial Hubble radius in conformal time. 

It is important to note that for the above approximation to work there must not exist 
fiow of modes from superhorizon to subhorizon regime, since for these modes the difference 
between the actual mode function //j, and the adiabatic mode function ^^f will not be small 
(recall that superhorizon modes evolve non-adiabatically) . Therefore, the above formulas 
can safely be used only for accelerating cosmologies. However, by comparing (26) with 
the known results, we will show below that there is good agreement even for decelerating 
cosmologies. In these cases the a parameter compensates the errors and the method works 
fine even for decelerating cosmologies. 

To check the validity of our method assume that the background evolution is fixed, i.e. 
the scale factor a{r]) is given. Then, (26) dictates how vacuum energy density evolves in 
time, which can be compared with the known results in the literature. 

Let us start with de-Sitter space and take 

a — — . (27) 

■n 

The vacuum energy density of a minimally coupled massless scalar field in the Bunch-Davies 
vacuum was calculated in [24], which can also be found in the book [1]: 

-29/^4 

= 960^- ^^^^ 

Integrating (26) we find 

(-119 + 60«^)/i^ Co 
960^^^^^ + ^^^^ 

Therefore, for Co = and a — (3/2)^/^ our calculation exactly matches the well-known 

result in de Sitter space. Since Bunch-Davies vacuum is de Sitter invariant the vacuum 
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— oo < n < 


Decelerating 


n = 


Minkowski 


< ra < 1 


Big-crunch 


n = 1 


de Sitter 


1 < n 


Accelerating 



TABLE I. Classification of the universes with power law expansion in the conformal time. The 
scale factor (30) in the proper time is given by o = (t/to)"^, where m = n/{n — 1). 

energy density must be a constant, which would justify the choice Co = 0. On the other 
hand, a parameter turns out to be of order one as expected. 
Consider now a background with power-law expansion 



(30) 



Using point-splitting regularization, the stress-energy tensor of a minimally coupled massless 
quantum scalar field in the Bunch-Davies vacuum has been calculated in [25], where the 

vacuum energy density is given by 



Pbd = 



28 807r2 
1 



D 



+ 



138247r2 



11527r2 
24nR + 24i?i?o° + 3i?^ 



\nl^]+ -0(2 + n)+ip{l-n) + ^ 



R 



Here is the di-gamma function, fx is an arbitrary constant and 



(3) if, 



-RR 



111/ ^^fi ^^pu 2 /^'^ ~ ~ ' 



QfivRpaR^'^ + -rR^g^u 



3 
(31) 

(32) 
(33) 



Scaling /i is equivalent to finite renormalization of a coupling constant in the quantum 
effective action (see e.g. comments below eq. (7.52) in [1]). Therefore, one may actually 
add an arbitrary multiple of ^^^i?o° to Pbd- In this way, the extra divergences that may 
arise when the argument of the di-gamma function is zero or negative integer can also be 

removed. 

Using (30), we evaluate (31) as 



Pbd = 



19207r%2a4 
90(n^ - 1) (in 



81 - 30n - 10971^ 

6{n + l)h^ 



np,'' 



-^-0(1 - n)+tl){n + 2) 



-1. 



(34) 
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When n — —1, which corresponds to a radiation dominated universe, ^^^Hq^ vanishes iden- 
tically, so there does not arise any issue with the logarithm and the di-gamma functions. In 
that case, the vacuum energy density can be found as 

= 96^' " = ^^'^ 
To avoid the subtleties related to the singularities of the di-gamma function, it is safe to 
take n to be a non-integer. 

The radiation dominated background n — —1 is special since the exact mode function can 
be determined from (7) as jjik — e"*'^^/ V2k, where the proper normalization is chosen for the 
Bunch-Davies vacuum. One then finds that (12) and (13) vanish identically. On the other 
hand, the finite adiabatic subtraction terms (14) can de determined as pp — — /i^/(9607r^a^), 
which gives 



geoTT^a^ 

The agreement between (35) and (36) shows the necessity of keeping finite adiabatic sub- 
traction terms for massless fields. In terms of our approximation for evaluating the integrals 

(12) and (13), eq. (36) corresponds a = 0, which may be considered as a singular case. 
Integrating (26) for the background (30) we find 



Pv 



2 



19207r2n2a4 
- 90(71^ - 1) In 



— (-21 + 153n + (79 + 60Q;^)n^ - (92 + 60q;^ + 60a^)n^) 



n-2 



^0 



By comparing (37) with (34), the Bunch-Davies vacuum can be seen to correspond Cq = 0. 
The coefficient of ln[h^) in (37) is a-indepcndcnt and this term agrees exactly with (34), 
which is highly non-trivial since one cannot adjust any free parameters. The other terms 
can actually be made to agree for any value of a by choosing the ratio /i/ho appropriately 
(as noted above this might be interpreted as a finite renormalization of a coupling constant 
in the effective action). 
One may insist to set 

11 = ho, (38) 

since both of these parameters are infrared cutoffs. In that case, to be able to get pv = Pbd 
the following equation must be satisfied: 
"36(n + 1)^ 
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+2^(l-„)+2^(„+2) = ^ 3(n^-l)(»-2) ' '^^^ 
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FIG. 1. The numerical plots of a-parameter as a function of the power n from (39). For larger 
\n\ the plots look very similar to (3,4) or (—4,-3) intervals for positive and negative values, 
respectively. 

It is not possible to solve a from this equation analytically. In Fig. 1 we give numerical plots 
for different intervals of unit size (as noted below (33), for integer values of n an infinite 
renormalization is required to make sense of pbd)- As n — ±oo, (39) reduces to 

6 ln(aV6) + (5 + 8a^ + 8a^) = 2^(1 - n) + 2ip{n + 2). (40) 

Using il){x) — '4){x — 1) = l/(x — 1), one sees that the solution for a in the interval (n, n + 1) 
becomes independent of n for large |?7.|. From the graphs it is obvious that a remains to 
be of order unity for generic values, which shows the validity of our approximation. On the 
other hand, for n G (1,2) and n G (—1,0) there turns out to be no solution for a. This 
last fact together with the observation that as n approaches to an integer a increases and 
in the limit ceases to be of order unity signal the need for a finite renormalization such that 

When n G (—1, 0), if one chooses /i = IO/iq, which corresponds to a — > 10a in (26), there 
exists solutions for a. Similarly, for n G (—1, 0), if one chooses /x = IOO/iq, which corresponds 
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FIG. 2. The numerical plots of a-parameter as a function of the power n for the interval (—1,0), 
where ^ = IO/iq and for the interval (1, 2), where /i = IOO/iq- 



to a — )■ 100a in (26), a becomes soluble. Corresponding numerical plots are given in Fig. 
2. Note that for any given n in these intervals there appears two solutions of a. 

In summary, we see that using our approximation it is possible to reproduce the known 
results in the literature. In de Sitter space (n = 1), it is enough to choose a = (3/2)^/"^. 
The radiation dominated universe (n = — 1) is singular since the result is given solely by the 
finite adiabatic subtraction terms, which corresponds to the choice a = 0. These are the only 
two cases with ^^^Hq'^ = 0. For all other integer values of n, one needs to apply an infinite 
renormalization to (34) by scaling the constant fi. We see that the logarithmic term, which 
is not affected by finite renormalizations, matches exactly in (34) and (37). The others with 
non-integer n can make to agree either by scaling fi or by adjusting a. It turned out that 
the approximation also works for negative values of n, i.e. for decelerating cosmologies, see 
Fig. 1. As discussed above, for decelerating backgrounds extra care is needed for the modes, 
which was born as superhorizon and becomes subhorizon in time. The error appeared for 
neglecting these modes is compensated by the a-parameter. 



IV. APPLICATION TO SLOW-ROLL INFLATION 

In this section we would like to apply our findings to the simplest inflationary scenario 
realized by a single infiaton field ip with a suitable potential V{(f). Our aim is to see 
whether the backreaction of a free massless quantum field has any impact on the standard 
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inflationary results. First, let us review the well-known results by ignoring the backreaction 
effects. As usual, the stress-energy tensor of the inflaton field is given by 

P^^\^' + V{^), (41) 
P^=\^'-V{v). (42) 

Recall that the dot denotes differentiation with respect to the proper time t defined as 
dt — adrj. The inflaton obeys the scalar equation 

ifi + SHip + d^V = 0. (43) 

To realize inflation in the slow-roll regime, one must assume that 

^^^V, < SHip. (44) 

To satisfy these conditions it is enough to take the slow-roll parameters and 77^ to be 
small: 

In this regime the inflaton equation becomes 
Moreover, the Friedmann equation implies 

If one defines the slow-roll parameters characterizing the change in the Hubble parameter 
and the inflaton as 

H ~ -eH^, (p^-{r]- e)Hip, (49) 

it is possible to show that 

6 = 6^, f] = f]^. (50) 

Note that e and 77 determine the tilt in the spectrum of the curvature perturbation and they 
are observationally important parameters of inflation. 
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Consider now the backreaction of the quantum field 0. Prom the differential equation 
(26) satisfied by pv, one can estimate that 



(51) 



Indeed, in the exact de Sitter space the vacuum energy density in the Bunch-Davies vacuum 
is given by (28). In the slow-roll regime one would not expect a large deviation from (28). 
In the following, we will treat e^, 77,^ and 



7 = 



M2 



(52) 



as small parameters and keep only the leading order terms. Since we do not assume any 
coupling between and (/?, the infiaton equation does not change: 



' 3H ■ 



(53) 



On the other hand, from the Friedmann equation 

Stt 



3M^ 



{PV+ P^): 



(54) 



one sees that ~ H^M^ and thus pv ^ ^ip^- Using (25), (14) and (15) in Einstein's 
equations, keeping only the leading order terms in the slow-roll parameters we get 
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(55) 



Prom (53) and the Priedmann equation (54), we find 

p 
An 



0^ = ^6,i/^ 



(56) 



which is also true in the absence of the quantum scalar field (i.e. the presence of pv in (54) 
does not alter the relation between infiaton kinetic energy and H'^ since e<^7 is assumed to 
be neghgible). Using this last equation in (55) leads 



119 a* ^ Sttpv 



7, 



(57) 



_3607r 67r ' 3H^ 

where e is defined in (49). Similarly, taking the time derivative of (53) and using the 
definition of rj in (49) , we find that 

77 = 77^. (58) 
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Therefore the backreaction of a quantum scalar field modifies the slow-roll parameter e, i.e. 
it is not simply determined by the infiaton potential. 

It is easy to see that the change in e is actually very small. For example, considering a 
chaotic inflationary scenario with the potential V = ^rn^ip^, it is known that the cosmologi- 
cally relevant perturbations exit the horizon when ip ~ 3Mp — 4Mp and a phcnomenologically 
viable value of the infiaton mass is given by m ~ lO^^Mp. One may then find that ~ 10~^ 
and 7 ~ 10~^^, thus the modification due to backreaction effects is tiny. 

While the backreaction effects are negligible during the near exponential expansion, they 
might be important for the subsequent cosmic evolution. As noted above, (26) shows that 
the vacuum energy density, which is driven by the terms in the right hand side, is of the 
order of H^. However, the homogeneous piece decreases like which is slower than 

in radiation or dust dominated universes. Therefore, one expects pv to become Hq 
during the exponential expansion, where Hq is the Hubble parameter of inflation, and then 
to decrease like 1/a^; thus one has pv ~ Hq/o? following inflation (this is indeed what is 
observed analytically in [15] in the simplified setting). On the other hand, the initial energy 
density of the radiation just after inflation can be estimated as HIM'^, so the radiation 
energy density decreases hke pr ~ HlM^/a^. Comparing these two, we see that when the 
universe expands Mp/Ho times after inflation, the energy density of the vacuum catches up 
the energy density of radiation and then it starts dominating the universe. In the above 
mentioned chaotic inflationary model Mp/Ho ~ 10^ and the equivalence corresponds to a 
very early cosmic epoch, for example, well before nucleosynthesis. The only loophole in the 
above conclusion is that (26) is rigorously valid in an accelerating background. However, as 
we saw in the previous section (26) gives correct results even for decelerating backgrounds by 
adjusting the a-parameter. Therefore, our findings indicate that the vacuum energy density 
created during inflation must be taken into account for succeeding cosmic evolution. 

Of course, in a realistic model the evolution of the vacuum energy density must be studied 
by considering different stages carefully. For example, during the reheating or preheating 
stage the coupling of our scalar to the oscillating infiaton field may change the evolution 
of vacuum energy density substantially. Thus, it might be possible to recover the current 
standard cosmological picture by such modifications. 
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V. CONCLUSIONS 



Determining how quantum fields affect the cosmological evolution, in the semiclassical 
approximation where gravity is treated classically, is an important problem. Quantum field 
theoretical effects can be crucial in resolving the initial big-bang singularity, during inflation 
or in explaining the presently observed acceleration of the universe. Compared to quan- 
tization in a fixed background, the backreaction problem is much more difficult to study 
since the mode functions of the field cannot be solved and thus their contributions to the 
stress-energy tensor cannot be determined explicitly. We manage to evaluate these contri- 
butions approximately by considering subhorizon and superhorizon modes separately. For 
a massless field, subhorizon modes with wavelengths smaller than the Hubble radius evolve 
adiabatically and their contributions cancel out by the adiabatic subtraction terms. On the 
other hand, the contributions of the superhorizon modes to the stress-energy tensor obey 
an equation of state. In that way, we are able to obtain an explicit equation for the vacuum 
energy density. To have a better approximation, we introduce a parameter of order unity 
to quantify subhorizon/super horizon division. By adjusting that parameter we show that 
our approximation applied to a field propagating in a fixed background exactly agrees with 
the known results in the literature. In our approach there arc important subtleties related 
to field being massless. One should be careful about the finite adiabatic subtraction terms 
which survive the massless limit. Moreover, there also exists infrared divergences which 
must be treated suitably. 

Wc apply our findings to a simple slow-roll inflationary scenario realized by a single 
scalar field. Wc observe that quantum fluctuations does not change the background evolu- 
tion appreciably during the near exponential expansion. They only induce a small change 
in one of the slow-roll parameters. However, vacuum energy density accumulated during 
the exponential expansion can be important at later stages, since it decreases like 
which is slower than radiation or dust. This decrease corresponds to an equation of state 
parameter ui = —1/3, therefore vacuum energy-density and vacuum pressure do not alter 
the acceleration equation for the scale factor of the universe following inflation. 

Encouraged by the fact that we can reproduce the known results in the literature by our 
approximation, the present work can be developed in different ways. It would be interesting 
to understand the situation with decelerating cosmologies better. For that, the contributions 
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of the modes, which were born as superhorizon and later on become subhorizon, must be 
estimated. In that way, one can sharpen the conclusion that the vacuum energy density 
created during inflation becomes important at later stages of cosmic evolution. It is also 
worth to the consider the cosmological backreaction of a massive field. In that case, there 
exists two scales in the problem, the mass of the field and the Hubble parameter, and 
subhorizon/superhorizon separation of the modes is more complicated than that of the 
massless field (for instance if the mass of the field is larger than the Hubble radius, all 
modes evolve adiabatically) . Moreover, for the contributions of the superhorizon modes to 
the stress-energy tensor, it seems difficult to write down an equation of state. Thus for 
a massive field one may attempt to use numerical analysis to solve the integro-differential 
equation system. 
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